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Abstract 

Gutzwiller's semiclassical trace formula for the density of states in a chaotic system 
diverges near bifurcations of periodic orbits, where it must be replaced with uniform 
approximations. It is well known that, when applying these approximations, complex 
predecessors of orbits created in the bifurcation ( "ghost orbits" ) can produce clear sig- 
natures in the semiclassical spectra. We demonstrate that these orbits themselves can 
undergo bifurcations, resulting in complex, non-generic bifurcation scenarios. We do so 
by studying an example taken from the Diamagnetic Kepler Problem. By application 
of normal form theory, we construct an analytic description of the complete bifurcation 
scenario, which is then used to calculate the pertinent uniform approximation. The 
ghost orbit bifurcation turns out to produce signatures in the semiclassical spectrum 
in much the same way as a bifurcation of real orbits would. 
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1. Introduction 



Since its discovery in the early 1970s, Gutzwiller's trace formula ^ has become a widely 
used tool for the interpretation of quantum spectra of systems whose classical counterpart 
exhibits chaotic behaviour. It represents the density of states of the quantum system as a 
sum over a smooth part and fluctuations from all periodic orbits of the classical system, 
which are all calculated from purely classical data. This formula assumes all periodic or- 
bits to be isolated, so that it fails close to a bifurcation of periodic orbits, where different 
orbits approach one another arbitrarily closely, leading to a divergence of the periodic orbit 
contributions. 

This failure can be overcome with the help of uniform approximations 2 ^' 5 ^ which take 
into account the contributions of all bifurcating orbits collectively. Although in generic 
Hamiltonian systems only codimension-one bifurcations can be observed as a single control 
parameter is varied, in practical applications of uniform approximations bifurcation scenar- 
ios of higher codimensions must be included: 7 ) If a periodic orbit successively undergoes 
several bifurcations, uniform approximations capable of collectively treating all orbits par- 
ticipating in the whole complicated bifurcation scenario are needed. 

In the construction of most uniform approximations known, ghost orbits 8 ) in the com- 
plexified classical phase space play a crucial role, because real periodic orbits that are born 
in a bifurcation tend to have ghost orbit predecessors before they turn real. So far, these 
ghost orbits have never been observed to undergo bifurcations themselves, although there is 
no a priori reason why this should be impossible. In fact, we will now present an example 
of a ghost orbit bifurcation which occurs in connection with a generic period-quadrupling 
bifurcation of a real orbit. We will show that the construction of a uniform approximation 
requires the inclusion of this bifurcation and that traditional normal form theory can be ex- 
tended to also cover bifurcations of ghost orbits. A more detailed presentation of our results 
can be found in 9), 10). 



§2. The bifurcation scenario 

As an example, we study the hydrogen atom in a magnetic field, which has been described 
in detail, e.g., in Refs. 11) -13). We assume the nucleus fixed and regard the electron as 
a structureless point charge. Due to a scaling property of the Hamiltonian, the classical 
dynamics does not depend on the energy E and the magnetic field strength 7 separately, but 
only on the scaled energy E = j^^E. To plot periodic orbits, we use scaled semiparabolical 
coordinates /i 2 = f + z, v 2 = r — z. 
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To look for ghost orbits, we complexify the classical phase space by allowing coordinates 
and momenta to assume complex values. Since our Hamiltonian is real, the system is sym- 
metric with respect to complex conjugation. Therefore, the complex conjugate of a periodic 
orbit is a periodic orbit itself, so that orbits usually occur in complex conjugate pairs. In 
exceptional cases, however, a periodic ghost orbit can coincide with its complex conjugate. 
These symmetric ghost orbits then have real periods and actions. 

One of the shortest periodic orbits of 
the diamagnetic Kepler problem is the bal- 
loon orbit, which we will now focus atten- 
tion on. It undergoes a period-quadrupling 
at a scaled energy of (in scaled atomic units, 
which we use throughout) E c = —0.342 025. 
If E > E c , a stable and an unstable real 
satellite orbit of quadruple period exist. At 
E c , they simultaneously collide with the bal- 
loon orbit and vanish. The real orbits are 
shown in figure 1 at the scaled energy of 
E = —0.340. The solid and dashed curves 
represent the stable and unstable satellite or- 
bits, respectively. For comparison, the bal- 
loon orbit is shown as a dotted curve. If 
E < E c , no real satellites are present, but there are a stable and an unstable ghost orbit 
instead, both of which are symmetric with respect to complex conjugation. Thus, at E c a 
generic island-chain bifurcation takes place. 

Furthermore, there is an additional ghost orbit in the phase space close to this bifurcation 
that is also symmetric with respect to complex conjugation. At E' c = —0.343 605 < E c , this 
orbit collides with the stable ghost satellite, the two orbits loose their symmetry and become 
a pair of complex conjugate ghosts at E < E' c . For scaled energies above and below E' c , 
the ghost orbits are depicted in figures 2 and 3, respectively. From the imaginary parts it 
can clearly be seen that the symmetry with respect to complex conjugation is lost in the 
bifurcation. 

This is the first example of a ghost orbit bifurcation that has been described in the 
literature. Its occurrence presents an additional challenge to the construction of a uniform 
approximation because an approximation that deals with the generic period quadrupling 
bifurcation only diverges at the bifurcation energy of the ghost orbit bifurcation and becomes 
undefined below E' c where the stable ghost satellite orbit used in its construction does not 
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Fig. 1. Real orbits at scaled energy E = 
—0.340 > E c , drawn in scaled semiparabol- 
ical coordinates. Solid and dashed curves: 
stable and unstable satellite orbits. Dotted 
curve: balloon orbit 
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Fig. 2. Real and imaginary parts of complex 
ghost orbits at scaled energy E = —0.343. 
Solid and dotted curves: stable and unsta- 
ble ghost satellites created in the period- 
quadrupling bifurcation of the balloon or- 
bit at E c = -0.342 025. Dashed curve: ad- 
ditional ghost orbit created in the ghost 
bifurcation at E' c = -0.343 605. 
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Fig. 3. Ghost orbits at scaled energy E = 
-0.344 < E' c . Solid and dashed curves: 
asymmeric ghost orbits (real parts coin- 
cide) created at the ghost bifurcation at 
E' c = -0.343 605. Dotted curve: unsta- 
ble ghost satellite created at the period- 
quadrupling bifurcation of the balloon or- 
bit 



exist any more. 

All periodic orbit parameters required for the construction of the uniform approximation 
were calculated numerically. The complete numerical data is described in detail in 9), 10). As 
an example, the orbital periods are shown in figure 4. From the successive confluences of the 
periods, the sequence of bifurcations becomes clearly visible: The period of four repetitions 
of the balloon orbit, which is always real, is indicated by a nearly horizontal line at T 5.84. 
Above E C) there are two additional solid curves representing the periods of the stable (upper 
curve) and unstable (lower curve) real satellite orbits. At E c , these curves change from solid 
to dashed as the satellite orbits become ghosts. Below E c , the unstable ghost satellite does 
not undergo any further bifurcations in the energy range shown, whereas the stable satellite 
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collides, at E' c , with the additional ghost orbit. The latter can clearly be seen not to be 
involved in the bifurcation at E c . Below E' c , these two orbits are complex conjugates of 
each other. Thus, the real parts of their periods coincide, whereas the imaginary parts are 
different from zero and have opposite signs. 
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Fig. 4. Orbital periods of the orbits involved in the bifurcations as functions of the scaled energy 
E = £ 7 - 2 / 3 



§3. The general form of the uniform approximation 

Before we return to classical normal form theory in section 4, we introduce, in this section, 
the basic formulas for the quantum density of states necessary for the construction of the 
uniform semiclassical approximation. The density of states of a quantum system with the 
Hamiltonian H can be expressed with the help of the Green's function G(E) = (E — H)~ l 
as 

d(E) = — — ImTr G(E) , (3-1) 

7T 

where the trace of the Green's function can be evaluated in the coordinate representation, 

Tr G(E) = J dx'dx 5{x' - x) G(x'x, E) . (3-2) 

The first step in the formulation of periodic orbit theory 1 ' 1 is to replace the Green's function 
G(x'x, E) with its semiclassical Van Vleck-Gutzwiller approximation. For systems with two 
degrees of freedom the semiclassical approximation to the Green's function reads 

G(x'x,E) = ]== V v^exp{^-S{x'x,E) - i-u) . (3-3) 

ihVlirih class . traj . U 2 J 

Here, the sum extends over all classical trajectories running from x to x' at energy E, S is 
the action of a trajectory, v its Maslov index, and D is given in terms of second derivatives 
of the action. The contribution of a single orbit to the trace can be evaluated by introducing 
coordinates parallel and perpendicular to the orbit. The integration along the orbit can 
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then be performed in a straightforward fashion. If one notes that on a periodic orbit the 
action S(x'x, E) is stationary with respect to the transverse coordinates, Gutzwiller's trace 
formula for isolated periodic orbits is finally obtained by integrating over the transverse 
coordinates using the stationary-phase approximation. It is this last step which fails close 
to a bifurcation, where periodic orbits are not isolated. 

The basic idea of the stationary-phase approximation is to approximate the action func- 
tion by a quadratic function in the neighbourhood of any individual periodic orbit. To 
achieve a collective treatment of the bifurcating orbits, we need to find an ansatz function $ 
on a Poincare surface of section which has got stationary points corresponding to all periodic 
orbits to be included in the uniform approximation. We can then relate $ to the classical ac- 
tion function S(x'x, E) by means of a suitable, albeit unknown, coordinate transformation. 
In terms of the uniform approximation can be shown to assume the form 3 ) 

= ex p{^o( £ ) - y} / ^^y^ p y) exp{^,Py)} • (34) 

The coefficient X can be evaluated at the stationary points of <P, where it assumes the value 

X , {m}T J\E^¥\ . (3-5) 
yJ\TrM-2\ 

Here, T and M denote the period and the monodromy matrix of the corresponding classical 
orbit, Hess^ the Hessian determinant of the normal form, and the notation {m} is meant 
to indicate that this factor does not occur at the satellite orbits. 

§4. Normal-form description of the bifurcation 

A systematic means to construct an ansatz function <P is provided by normal form the- 
ory. 14 )' 15 ) Here, we adopt a normal form used by Schomerus 16 - ) to describe codimension-two 
bifurcations. In canonical polar coordinates (I,<p), which are connected to Cartesian coor- 
dinates (p, q) by 

p = V2I cos <p , q — V2I sin ip , (4-1) 

it reads 

<P = eI + aI 2 + bl 2 cos(4v?) + c/ 3 (l + cos(4y9)) . (4-2) 

This normal form turns out to qualitatively describe the sequence of bifurcations encountered 
here for suitably chosen parameter values. 

To establish the connection to the classical bifurcation scenario, we have to determine 
the stationary points of <P. The central periodic orbit at I — does not show up as a 
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stationary point because the polar coordiante system (4-1) is singular there. In the case 
\a\ > \b\ and c < we find three stationary points, each of which appears at four different 
angles ip because it corresponds to an orbit of quadruple period. For these stationary points, 
the real parts of the action coordinates / are shown in figure 5. The dotted line in this figure 
corresponds to a pair of complex conjugate stationary points. 

To interprete these results, we observe that according to its definition (4-1) the coordinate 
/ is positive for real orbits and that the action <P(I, ip) is real for real I, ip. Therefore, negative 
real solutions / correspond to ghost orbits which are symmetric with respect to complex 
conjugation and thus have real actions, whereas a complex / indicates an asymmetric ghost 
orbit. 

If e > 0, we have two stationary points 
at positive values of / and one stationary 
point at a negative I. They correspond to 
two real satellite orbits of quadruple period 
and a symmetric ghost orbit. As e decreases 
through zero, two stationary points simul- 
taneously move from positive to negative /, 
thus indicating that the two real satellites 
become ghosts in an island-chain bifurca- 
tion. Finally, two stationary points with 
negative action coordinates collide at a neg- 
ative value of e and become complex. This 
describes a ghost orbit bifurcation in which 
two symmetric ghost orbits collide and loose 
their symmetry. 

The normal form (4-2) thus qualita- 
tively reproduces the bifurcation scenario 
described above. The normal form parameters e, a, b, c have to be determined so as to 
make the description quantitatively correct. To this end we calculate the stationary values 
of the normal form, equate them to the actions of the periodic orbits and then solve for the 
normal form parameters. As can be shown by a lengthy calculation, 10 ) these parameters 
are uniquely determined by the actions and the energy if we choose e = E — E c and require 
all parameters to depend continuously on the energy. 
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Fig. 5. A sketch of the bifurcation scenario 
given by the normal form (4-2) for the case 
\a\ > \b\ and c < 0. Solid curves: real or- 
bits. Dashed curves: ghost orbits symme- 
tric with respect to complex conjugation. 
Dotted curve: a pair of complex conjugate 
ghosts. 
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§5. Evaluation of the uniform approximation 

After the ansatz function has been completely specified, a suitable approximation to 
the coefficient X in (34) remains to be found. We assume X to be independent of if, and as 
the value of X is known at the stationary points of <P at four different values of I (including 
I = 0), we approximate X by the third order polynomial p(I) interpolating between the four 
given points. This choice ensures that our approximation reproduces Gutzwiller's isolated- 
orbits formula if, sufficiently far away from the bifurcations, we evaluate the integral in 
stationary-phase-approximation. Thus, the uniform approximation takes its final form 

d{E) = ^Re exp{iSb(£) - i\i>] J dYdP{,p{I) exp {^>(F, P^)} , (5-1) 
which contains known functions only and can be evaluated numerically. 
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Fig. 6. Uniform approximation to the contribution of the considered bifurcations to the density of 



states for two different values of the magnetic field strength: (0)7 = 10 12 , (6)7 = 10 14 . Solid 
curves: uniform approximations. Dashed curves: Gutzwiller's trace formula. 

We calculated the uniform approximation (5-1) for two different values of the magnetic 
field strength 7. The results are shown in figure 6. To suppress the highly oscillatory 
contributions originating from the factor exp {jrS'o(P)}, we plot the absolute value of (5-1) 
instead of the real part. As can be seen, the uniform approximation proposed is finite at 
the bifurcation energies, and, as the distance from the bifurcations increases, asymptotically 
goes over into the results of Gutzwiller's trace formula. Even the complicated oscillatory 
structures in the density of states caused by interferences between the contributions from 
the different real orbits involved at E > E c = —0.342 025 are perfectly reproduced by our 
uniform approximation. We also see that the higher the magnetic field strength, the farther 
away from the bifurcation the asymptotic (Gutzwiller) behaviour is acquired. The magnetic 
field dependence of the transition into the asymptotic regime can be traced back to the 
fact that, due to the scaling properties of our system, the scaling parameter 7 1 / 3 plays the 
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role of an effective Planck's constant, therefore the lower 7 becomes, the more accurate the 
semiclassical approximation will be. 

§6. Conclusion 

We have shown that in Hamiltonian systems with mixed regular-chaotic dynamics ghost 
orbit bifurcations can occur besides the bifurcations of real orbits. These are of special 
importance when they appear in the vicinity of bifurcations of real orbits, since they turn 
out to produce signatures in the semiclassical spectra much the same as those of the real 
orbits. Consequently, the traditional theory of uniform approximations for bifurcations of 
real orbits must be extended to also include the effects of bifurcating ghost orbits. 

We have illustrated the phenomenon of bifurcating ghost orbits in the neighbourhood of 
bifurcations of real orbits by way of example for the period-quadrupling of the balloon orbit 
in the diamagnetic Kepler problem, and have demonstrated how normal form theory can 
be extended for this case so as to allow for a unified description of both real and complex 
bifurcations. 

We picked the example mainly for its simplicity, since (a) the real orbit considered is one 
of the shortest fundamental periodic orbits in the diamagnetic Kepler problem and (b) the 
period-quadrupling is the lowest period-m-tupling possible (m = 4) that exhibits the island- 
chain bifurcation typical of all higher m. Thus we expect ghost orbit bifurcations to appear 
also for longer-period orbits, and, in particular, in the vicinity of all higher period-m-tupling 
bifurcations of real orbits. 

In fact, a general discussion 10 ) of the bifurcation scenarios described by the normal form 
(4-2) and a more general variant of its for different values of the parameters leads us to the 
conclusion that the appearance of ghost orbit bifurcations in the vicinity of bifurcating real 
orbits is the rule, rather than the exception, in general systems with mixed regular- chaotic 
systems, and thus one of their generic features. 
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